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. To measure the strong coupling as from event shape observables two ingredients are 
necessary. A perturbative prediction containing the dependence of observables on as and 
>■ I a description of the hadronisation process to match the perturbative prediction with the 
^ I hadronic data. 

Pn| ' As perturbative prediction 0{a1), NLLA and combined calculations are available. Be- 
side the well known Monte-Carlo based models also analytical predictions, so called power 
Q>^ ! corrections, exist to describe the hadronisation. Advantages and disadvantages of the dif- 
G^ I ferent resulting methods for determining the strong coupling and its energy dependence 
will be discussed, the newest Delphi results will be presented, and an overview of the 
LEP results will be included. 

^ . 1. INTRODUCTION 



The theoretical description of hadron production in e~^e -annihilation consists of three 
parts. The first part is based on the fundamentals of the standard model: Feynman dia- 
grams are used to (perturbatively) calculate the electroweak process of e~^e -annihilation 
and the evolution of partons under the strong interaction. At some stage the partons 
must be combined to become hadrons. This hadronisation process cannot be described 
by perturbation theory and thus builds a second (non-perturbative) part. Finally the 
decay of unstable hadrons, which can be described by kinematics using experimentally 
measured decay rates, need to be included before the prediction can be confronted with 
data. 

Several different predictions exist for the two first parts. Calculations of the evolution 
of quarks and gluons are available in fixed order 0{a'^), recently some observables be- 
came available in 0{al), and in the next to leading log approximation (NLLA), which 
resums large logarithms to all orders of a^. When combining 0{al) with NLLA matching 
ambiguities occur leading to even more competing predictions. 

To describe the second part (the hadronisation) usually generator based models are 
used. The most reliable of these models are the Lund string fragmentation, implemented 
in Jetset, and the Cluster fragmentation, implemented in Herwig. Recently analytical 
predictions, so called power corrections, became available to describe the influence of 
hadronisation on event shape observables. 

Given the many different possibilities to determine it was necessary to restrict this 
review due to time and space limitation. I've chosen to review the analyses using event 
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shape observables to determine a^, as this field currently is very active at LEP. 

In the following I will first discuss the standard method for measuring from event 
shapes, which uses 0{a1), NLLA or the combined (9(a^)+NLLA prediction in combina- 
tion with generator based hadronisation models. In the second part the energy dependence 
of will be discussed introducing the alternative analytic description of hadronisation 
by power corrections. Finally first results obtained using the new C(a^) predictions are 
presented. In all three sections DELPHl-analyses shall serve as showcases. 



2. CONSISTENT DETERMINATION OF THE STRONG COUPLING 

The original motivation of repeating an a^-measurement with LEPl data in Delphi 
was to include the event orientation, given e.g. by the polar angle of the Thrust axis 
9ti into the fit. Second order coefficients including such an event orientation became 
available in 0{a1) through the Event2 program by Catani and Seymour Using the 
fully reprocessed LEPl data, event shapes were measured very precisely in eight bins 
of 9t- The small systematic errors result not only from the good quality of the final data 
reprocessing, but also from the fact that all detector corrections in this measurement were 
naturally calculated for each ^-r-bin separately. 

The second order prediction for such distributions now contains second order coefficients 
A and B which depend on the observables value itself and on 9t'- 

m - ^<-^^'^ ^ (^'-^^' ■ ^-*»'" (£) ^^'-^^ « 

/X being the renormalisation scale and 69 = (33 — 2iVj)/127r. Beside this formula 
contains renormalisation scale parametrised as = 11^ / E^^ as a free parameter. 

As a first step the dependence of the resulting a^-values on this parameter were inves- 
tigated. It is found that different observables aquire largely different dependencies on x^. 
Also the scales with the optimal x^/ndf vary widely (see Fig. ^left and middle). In spite 
of this large variation of the optimal scale, the results corresponding to the experimen- 
tally optimised scales show a smaller spread for a large number of different event shapes, 
than the results obtained with x^ = 1. In contrast to the results with scale 1, optimised 
scales yield consistent values of Ug without assuming any renormalisation scale error (see 
Fig. right). In addition the scale dependence of near the optimised scale is smaller 
than for a;^ = 1, so that the scale variation between half and twice the chosen value of x^ 
leads to a smaller scale uncertainty for optimised scales. 

Averaging over all 18 investigated observables the final result is 

a,(Mz) = 0.1228 ±0.0119 (a;^ = 1) (2) 
as{Mz) = 0.1173 ±0.0026 (opt. scales), (3) 

where for the optimised scales also the 6-quark mass corrections are included. The overall 
fit quality of these results is far better for the optimised scales than for x^ — 1, as can 
be seen from the Ax^ in Fig. |l]. Thus the larger spread of the results with = 1 can be 
attributed to a less stable fit procedure which is caused by the bad agreement between 
data and the prediction. 
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Figure 1. Dependence of the resulting as on the renormahsation scale for different event 
shapes (upper parts in left and middle plot) and the corresponding A^^-values (lower 
parts). Resulting as values for 18 different event shape observables using experimentally 
optimised scale and including corrections for fe-quark mass (right plot). 



To crosscheck the results the scale dependence was also investigated for NLLA and 
combined 0{al)+NLLA predictions. In contrast to the 0{a1) results the required relative 
renormahsation scales are close to one and thus there is no significant change compared to 
the results obtained with = 1. Moreover, because of the limited fit range available for 
NLLA and because of the matching ambiguities for combined (9(a;^)+NLLA predictions 
the total error for these methods is larger than for the O(a^) result: 

as{Mz) = 0.116 ±0.006 NLLA (4) 
as{Mz) = 0.119 ±0.005 0(a2)+NLLA (5) 

Both results are in good agreement with each other and with the 0{a'^) results 0. 

There are two other publication in which optimised scales are used to determine a<j at 
Ecm = Mz'- Already in 1992 Opal |Q states to have found a clear improvement in the fit 
quality with optimised scales using 14 observables. They quote 

a,(Mz) = 0.1181°:°°^ (opt. scales), (6) 

where the error includes a scale variation from the optimised scale upto a scale of 1. 

In 1996 Burrows et. al. investigated 15 Observables using SLD-data. In contrast to 
Opal and Delphi they found no significant reduction of the spread of a^-values, though 
the shift to lower values of when using optimised scales is reproduced: 

a,(Mz) = 0.1265 ±0.0076 (a;^ = 1) (7) 
a^(Mz) = 0.1173 ±0.0071 (opt. scales). (8) 
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In spite of the extra theoretical uncertainties due to matching ambiguities three of four 
LEP experiments today use the combined (9(a^)+NLLA calculation to determine their 
central a^-value (second error at L3 is the theoretical component): 

asiMz) = 0.1216 ±0.0039 Aleph [5] 

as{Mz) = 0.1220 ±0.0015 ±0.0060 L3 § (9) 
a,(Mz) = 0.120 ±0.006 Opal 0. 

3. ENERGY DEPENDENCE 

The increase of beam energy accomplished during the LEP2 programme gives access to 
the energy dependence of event shapes and thereby to the energy dependence of a^. 

3.1. Power Corrections 

Using data at different energies allows also to replace the generator based hadronisation 
models by an analytical ansatz. For mean values this ansatz (which was developed by 
Dokshitzer and Webber |^) describes the hadronisation by an additive term: 

(/) = —/ f^da = (/pert) ± (/pow) . (10) 
(Ttot J da 

The 2nd order perturbative prediction is given by Eq. (|l]) with coefficients A and B 
integrated over y and 6t- The power correction term is falling off like the inverse centre- 
of-mass energy and is given by 



ifpov.) = Cf—^M — 
-C/rm 



ao(/i/) - «s(/i) - ( 6o ■ log ^ ± ^ ± 26o j al{ii) 




(11) 



where ao is a non-perturbative parameter accounting for the contributions to the event 
shape below an infrared matching scale K = (67/18 — 7!'^/6)Ca — 5Nf/9. The 
Milan factor Ai is set to 1.8, which corresponds to three active flavours in the non- 
perturbative region ||^. The observable-dependent coefficient is 2 and 1 for / = (1 — T) 
and / = (M^/E'^jg), respectively. For (-Bmax) the coefficient is itself energy dependent: 
C(Bmax> ~ 1/ \J a s{E cm) 0- The infrared matching scale is set to 2 GeV as suggested by 
the authors 0, the renormalisation scale /i is set to be equal to Ecm- Beside these 



Observable 


ao(2GeV) 


as(Mz) 


xVndf 


(-^max) 

(1-T) 
{Ml/El,:^ 


0.407 ±0.022 ±0.055 
0.493 ± 0.009 ± 0.006 
0.545 ±0.023 ±0.017 


0.117 ±0.0012 ±0.0015 
0.119 ±0.0014 ±0.0067 
0.120 ±0.0020 ±0.0050 


8.28/21 
64.0/32 
8.28/21 



Table 1 

Determination of in a combined fit of ctg and using a large number of measure- 
ments [|10| at different energies |jl2|. The first error is calculated from the fit, the second 
contains the renormalisation scale dependence. 



^ 0.13 
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[£0.125 
^ 0.12 
0.115 
0.11 
0.105 
0.1 
0.095 
0.09 
0.085 



DELPHI prelim. 




0(a;) (Power Corr-.) 

Stalistical En"or.s only 

World average (PDG 1998) 
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O(a^) 

Statistiual Errors only 

World average (PDG 1998) 
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Pure NLLA 

Statistical Errors only 

World average (PDG 1998) 
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NLLA+0(a;) log(R)-Scheme 

Statistical Enors only 

World average (PDG 1998) 
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Power Corr. Monte Carlo based hadronisation models 

0(^2) NLLA C»(a2)+NLLA 



dlog£;c 



1.07±0.19±0.20 L23±0.31 ±0.19 



L41±0.31 ±0.27 



1.12±0.22±0.18 



Figure 2. Comparison of a^-results obtained by Delphi at different centre-of-mass ener- 
gies using Monte Carlo generator based hadronisation models as well as power correction. 
Below each plot the result for the logarithmic energy slope of the inverse coupling is given. 



formulae contain ao as the only free parameter. In order to measure as from individual 
high energy data this parameter has to be known. 

To infer ao, a combined fit of as and ao to a large set of measurements at different 
energies fl^ is performed |]lll|l2|. For Ecm > Mz only Delphi measurements are included 



in the fit. The resulting values of ao for (-Bmax), (1 — T) and (Ml/E^^^) are summarised 
in Tab. |I[ Even though the found ao values are experimentally inconsistent, the univer- 
sality of ao is not violated because the expected theoretical precision allows deviations of 
upto 20%. The a^ values are consistent with the corresponding analyses of L3 [|] and 
Opal/Jade |T3[ . 



After fixing ao to the values found, a^ can be calculated individually for each energy 
using Eqs. (p!0|-[Tl|). The results for energies between 65 GeV and 189 GeV of this method 
and of the traditional methods described in the previous section are compared in Fig. ^. 
All methods give consistent results. 

3.2. Energy dependence of ocg 

To measure the energy dependence of a^ DELPHI uses the logarithmic energy slope 
of the inverse coupling. This quantity is directly proportional to the Callan-Symanzik 
/3-function and is independent of a^ and of Ecm to first order: 

da-^ 1 Po (3_ 



-^/5«. = ^ + fV«^ + .-.^l-27 (12) 



dlogEem a^ 27r An^ 



The numerical value represents the QCD prediction calculated in second order for energies 
between 91 GeV and 200 GeV using the PDGs world average of ag. The energy dependence 
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of this derivative in the given range and the uncertainty of as influence this value by about 
one unit in the last digit. 

The result with the smallest systematic error is obtained from 0{a^^)+NLLA fits: 

= 1.12 ±0.22 ±0.18 (13) 

d log Ecm 

in good agreement with the QCD expectation. 

Instead of determining a value for the energy dependence explicitly L3 checks the 
running of by applying a combined fit to all energies assuming the standard model 
running to C(a^). This yields a x^/ndf = 15.4/12 also indicating a good agreement 



4. 0{al) 

Third order calculations for four-parton final states have become available recently [|n 
[15|. These calulations can be used to measure in next to leading order from event 
shapes that acquire non-trivial values only for four and more partons, like the four-jet- 
rate R^. 



cm 



Observables of this kind are uncorrelated to the observables discussed so far which are 
based on three-jet-like configurations. The reduced number of relavant events is partly 
compensated by the (due to quadratic dependence) increased sensitivity in Eq. (|1^). 

Delphi investigated four-jets-rates RA^ycut) at a given ycut for the Durham cluster 
algorithm |l^. It shows good agreement with the prediction for ?/cut > 0.002. Using 
Ucut = 0.0025 the resulting as values are consistent with the results shown in the previous 
section, but they have larger statistical errors. Also the running obtained from this 
analysis is in good agreement: 

da7^ , , 

= 1.16 ±0.46 . (15) 



dlogE, 



This new calculation has also been investigated by Aleph, but was not yet used to 
determine the strong coupling. 



5. SUMMARY 

In this talk it was tried to give an overview over the current state of the art in measuring 
as from event shape observables at LEP. Beside different perturbative calculations several 
hadronisation models (generator based and analytical ones) exist. The previous standard 
choice for the perturbative calculation 0(q;^)±NLLA is questioned by Delphi due to 
poor data description and unsatisfactory consistency. 

In that sense measurements of from LEP data are still in full progress. In addition 
to the method used, special attention is payed to tests of the running, also new theoretical 
developements like power corrections and 0{p?^ gain due recognition. 
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